Collision geometry scaling of Au + Au pseudorapidity density 
from = 19.6 to 200 GeV 
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The centrality dependence of the midrapidity charged particle multiplicity in Au + Au collisions 
at y/SNN ~ 19.6 and 200 GeV is presented. Within a simple model, the fraction of hard (scaling 
with number of binary collisions) to soft (scaling with number of participant pairs) interactions 
is consistent with a value of a; = 0.13±0.01(stat)±0.05(syst) at both energies. The experimental 
results at both energies, scaled by inelastic p{p) + p collision data, agree within systematic errors. 
The ratio of the data was found not to depend on centrality over the studied range and yields a 
simple linear scale factor of R200/19.6 = 2.03±0.02(stat)±0.05(syst). 

PACS numbers: 25.75.-q, 25.75.Nq 



We have studied the centrality dependence of the 
charged particle multiplicity at midrapidity for Au + Au 
collisions at nucleon-nucleon center of mass energies 
^/sWn = 19.6 GeV and 200 GeV, using the PHOBOS 
detector. The data were taken during the Au-\-Au run in 
2001-2002 at the Relativistic Heavy Ion Collider (RHIC) 
in Brookhaven National Laboratory. Data at both ener- 
gies have allowed the extraction of results with the same 
detector, which covers a factor of ten in collision energy, 
from slightly above the highest energy of the CERN SPS 
fixed target program to the highest RHIC energy. 

The latest results from RHIC have suggested the ef- 
fect of 'jet quenching' in central Au -\- Au collisions that 
acts to reduce both the overall yield of high pr parti- 
cles Q, H, il and back- to-back jet correlations The 
presence of these dramatic effects for the most central 
Au -\- Au collisions at y^s^N = 130 and 200 GeV, as 
well as their absence in the cases of peripheral Au -\- Au 
0, H , central d + Au and inclusive d Am HQ 

have been generally well reproduced by calculations that 
utilize a pQCD framework to calculate the initial high p^ 
production rates, coupled with a large energy loss in the 
dense medium |9|. In this picture, the produced 'dense 
medium' is responsible for the experimental effect, which 
presumably occurs only in the large overlap volume of 



central Au -\- Au collisions. 

One of the intriguing overall features of the RHIC data, 
however, is that models solely based on parton satura- 
tion in the colliding nuclei describe the detailed centrality 
and rapidity dependence of the measured charged parti- 
cle multiphcities at 130 and 200 GeV [ll[I3- If parton 
saturation is playing a significant role in these relativis- 
tic heavy ion collisions, it could also reduce the initial 
production rate of high p^ particles to the extent that a 
large energy loss in the dense medium is no longer nec- 
essary to describe the data for central Au -\- Au collisions 
0, . Since the measured total charged particle multi- 
plicities are completely dominated by the emission of low 
Pt (< 1-5 GeV/c) particles, one might speculate that the 
production dynamics of low p-r particles is quite different 
from those at high p^- A study of the detailed centrality 
dependence of the bulk charged particle production over 
a large energy range may, therefore, provide additional 
constraints on models attempting to describe both the 
low and high pT behavior of particle production. 

The PHOBOS detector configuration was the same for 
measurements at ^/SNN — 19.6 and 200 GeV. Specifically, 
the detectors used in this analysis were the centrally lo- 
cated Octagon barrel. Vertex detector and the multiplane 
Spectrometer. These detectors are all constructed from 
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FIG. 1: Pseudorapidity density distributions from y^ijvjv ~ 
200 (light band) and 19.6 (dark band) GeV Au+Au collisions, 
for the most central 25% of the cross section 171]. The boxed 
areas (a-d) illustrate the separate regions in pseudorapidity 
used for the centrality determinations at each energy. 



silicon wafers, more details can be found in Refs. [l^ll3| . 
The primary trigger for the data reported here is based 
on 71 > 2 hits in two segmented, large-area scintillator 
counter arrays (Paddles) covering 3.2 < jr/l < 4.5 rela- 
tive to the nominal vertex position. Pseudorapidity is 
defined as 77 = — lntan(6'/2) where 9 is the polar angle 
to the beam axis. This trigger was sensitive to 88% of 
the total inelastic cross section in the 200 GeV data [l^ 
and is estimated to be 80% efficient for 19.6 GeV, from 
Monte Garlo (MC) studies using Hijing JJi] and a full 
GEANT [3 simulation of the PHOBOS detector. 

The centrality determination for the PHOBOS Au+Au 
results at = 130 and 200 GeV [HEB uses the en- 

ergy signals from the Paddle counters, which lie away 
from midrapidity, as illustrated by region (b) in Fig. ^ 
These signals were found in the MC simulations to in- 
crease monotonically with the number of nucleons par- 
ticipating in the collision and therefore, through bins in 
the percentage of total cross section, provide a measure 
of centrality. A similar centrality measure can be cre- 
ated at 19.6 GeV by scaling the Paddle pseudorapidity 
range by the ratio of beam rapidities at 200 and 19.6 
GeV, yfgl^/yfoo"" = 0.563, as shown in region (d) of 
Fig. HI The resulting 7] region, 1.8 < |?7| < 2.5, lies within 
the Octagon detector coverage of I77I < 3.2 for coUisions 
which occur within ±10 cm of the nominal vertex posi- 
tion. Thus, a centrality measure based on the deposited 
energy in this specific region (d) of the Octagon was cal- 
culated for 19.6 GeV and used in determining centrality 
for direct comparison to the original Paddle based cen- 
trality determination at 200 GeV. 



A second, independent, centrality measure was also de- 
veloped for both the 200 and 19.6 GeV datasets. This 
was based on the energy of charged particles traversing 
the Octagon within |?7| < 3.0 at 200 GeV, region (a) of 
Fig. n and a 'reduced' region of \t]\ < 1.8 at 19.6 GeV, 
region (c). 

The vertex position of each event is required for the 
merging and angle correction of valid hits in the Oc- 
tagon. The primary collision vertex used was determined 
by straight-line tracks in the first six planes of the Spec- 
trometer. The same vertexing algorithm was used at 
both energies. From MC studies of the detector, this 
vertex has a resolution (Jx.y,z ~ 0.3,0.3,0.4 mm for cen- 
tral and <7x,y,z ~ 0.6,0.5,0.8 mm for mid-peripheral col- 
lisions at 19.6 GeV, where z is along the beam and y 
is vertical. For more peripheral collisions, the efficiency 
falls away smoothly from 100% with decreasing central- 
ity. Additional cross-checks performed with different ver- 
texing methods yielded consistent results. 

Due to the requirement of a valid vertex, the resulting 
dataset is not only biased by the intrinsic trigger effi- 
ciency, but also by our (track-based) vertex reconstruc- 
tion efficiency. The vertex biased detection efficiency at 
200 GeV is deduced from prior measurement 
19.6 GeV, the efficiency is determined from a ratio of 
yields (Data/MC) after shape matching of multiplicity 
distributions between data and MC simulations in the 
centrality regions (c) or (d) from Fig. ^ The match- 
ing algorithm only utilizes data where 100% efficiency in 
both triggering and vertexing is expected. The overall 
efficiency, including the trigger and vertexing bias, is es- 
timated to be 66.0 ± 2.0 % and 55.4 ± 2.0 % for the 200 
and 19.6 GeV data, respectively. 

The final step in the centrality determination is to con- 
nect the experimentally deduced cross section percentiles 
with a well-defined variable, such as the number of partic- 
ipating nucleons, Npart- In order to do this, a monotonic 
relation was assumed to exist between the multiplicity 
distribution in the chosen 77 region and Npart- This as- 
sumption has been borne out by extensive simulations 
and the experimental (inverse) correlation between mul- 
tiplicity and zero degree calorimetry. Additionally, in or- 
der to further reduce ambiguity in this assumption, the 
measures of centrality outlined above were chosen specif- 
ically to lie in very different regions of pseudorapidity. 
Away from midrapidity, regions (b) and (d), particle pro- 
duction in the Hijing model depends linearly on Npart, as 
discussed in Ref. |l8l |. However, the midrapidity charged 
particle production, regions (a) and (c), is not linear as 
predicted by Hijing. Comparison of the results from these 
distinct regions for centrality determination can expose 
any systematic effects of non-linearity. 

The 'tracklet' reconstruction method was used to ob- 
tain the best measure of the charged particle multiplic- 
ity at midrapidity. A tracklet is a two-hit combination 
from the inner and outer Vertex detector layers, which 
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FIG. 2: The measured pseudorapidity density per participant 
pair as a function of (Npart) for Au + Au collisions at ^ysNN ~ 
19.6 GeV (closed squares), 200 GeV (closed circles). The 
error ellipses around the data combine the systematic error 
on dNch/dril^^^^i and (Npart)- The open circle and square 
represent the p{p) +p results (see text). The three curves 
give two model calculations and one fit result. As in [in| |. the 
systematic errors represent 90% C.L. limits. 



points back to the reconstructed vertex. A Vertex detec- 
tor tracklet is only formed when the difference between 
the hits (residuals) in azimuthal angle and rj are less than 
0.3 radians and 0.04 respectively. The difference in the 
magnitude of these values originates from the granular- 
ity of the detector in the respective measured directions. 
The final multiplicity is measured in the region jryj < 1 
by counting all reconstructed tracklets and correcting 
for combinatorial background, detector acceptance and 
background from secondary particles and weak decays. 
The correction factors, determined using MC simulations 
as a function of reconstructed vertex position and num- 
ber of hits in the vertex detector, were found to be the 
same for each energy to within 1%. The results of the 
Vertex detector tracklet analyses obtained using central- 
ity regions (a) and (c) of Fig. are given in Table n We 
obtain the same result for the yield per participant pair, 
dNch /di]/{^ (Npart)), from centrality regions (b) and (d). 
Systematic errors for the Vertex detector tracklet results, 
as in [13, are 7.5% (90% C.L.) at both energies. 

An independent analysis was also carried out at both 
energies using the tracklet method applied to data ob- 
tained from pairs of hits in adjacent Spectrometer planes. 
This additional analysis yielded results consistent within 
2% across all centralities to the results obtained using 
the Vertex detector. 

The data for 19.6 GeV, together with the reanalyzed 
results for 200 GeV, are shown in Fig. |21 This new 
result for 200 GeV has a slightly flatter dependence 



on (Npart) than found previously |lO|, but within the 
quoted systematic errors. This flattening of the yield 
with (Npart) for central collisions arises entirely from the 
new (vertex restricted) centrality measures and methods 
detailed above. Also shown, at Npart — 2, is the in- 
elastic charged particle multiplicity obtained in p{p) + p 
collisions ^^ 0| for data measured at 200 GeV 
{dNch/dr]\\r,\<i = 2.29 ±0.08) and interpolated for 19.6 
GeV idNch/dr]\\rj\<i = 1.27±0.13). The p(p) -hp data are 
averaged over the same pseudorapidity region, \r]\ < 1, as 
the heavy ion measurement. Clearly, the yield of charged 
particles per participant pair at midrapidity for the mea- 
sured Au + Au collisions is higher than found in corre- 
sponding p{p) -\- p collisions. 

The dotted line in Fig. [3 represents a fit to the data 
using the sim ple two-component parameterization pro- 
posed in Ref. 
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The value for Ncou, the number of binary (nuclcon- 
nucleon) collisions, is determined from a Glauber model 
calculation and is found to depend on the number of par- 
ticipants, Npart, through a simple power law. We find 
that Ncoii = Ax N^art, with A = 0.33 and 0.37 and a 
= 1.37 and 1.32 for 200 GeV and 19.6 GeV respectively. 
The difference in the A and a parameters at the two 
energies is due to the different measured nucleon-nucleon 
cross sections at y^J/vw = 200 GeV (aNN = 42±1 mb) vs. 
19.6 GeV (ajsiN = 33±1 mb). The remaining parameters 
are Upp, the yield obtained in p{p) + p collisions, and x, 
which represents the contribution from 'hard' processes 
taken to scale with Ncoii- 

The large systematic errors on the data preclude a si- 
multaneous extraction of both Upp and x solely from the 
Au + Au data. If only statistical errors are considered, 
the extracted parameters are npp « 2.7, 1.3 and x ~ 0.09, 
0.11 at 200 and 19.6 GeV, respectively. 

A value for x can also be obtained by fixing Upp at 
the measured and interpolated values of 2.29 and 1.27. 
Using statistical errors, we find x = 0.145 and 0.120 for 
200 and 19.6 GeV, respectively, as depicted by the dotted 
line in Fig. [21 The systematic error on the fit parameter 
X was determined by allowing the pp value and the data 
points to vary independently within their systematic un- 
certainties. Within the systematics, we find the fraction 
of hard collisions for both energies is consistent with a 
single value of x = 0.13±0.01(stat)±0.05(syst). 

The equivalence of parameter x at both energies, 
within the large errors, is surprising as the pQCD cross 
section for processes with large momentum transfers is 
expected to rise from ^/snn = 19.6 to 200 GeV. This ex- 
pectation of increasing slope in centrality with collision 
energy is shown by the dashed lines in Fig. El which rep- 
resent Hijing predictions [l^ . Although this rapid rise 
is not conclusively ruled out within the systematics at 
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TABLE I: Experimental results for the charged particle pseudorapidity density at midrapidity as a function of percentile 
cross section. The most central collisions are labeled as Bin 0-3%. The requirement of 100% efficiency in both triggering and 
vertexing imposes a lower limit on the reported results of (Npart) > 65 and 95 for 200 and 19.6 GeV, respectively. All errors 
represent 90% C.L. systematic limits with the exception of the Ratio, for which the errors are standard combined statistical 
and systematic 1-a uncertainties. 
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y^SNN = 200 GeV, it clearly does not follow the trend 
found in the data. Calculations from the parton satura- 
tion model 0, [2^ (solid lines in Fig. El predict a much 
weaker centrality dependence for both energies, in better 
agreement with the experimental data. 

In order to gain a different perspective on the data, we 
scale the Au -\- Au charged particle pseudorapidity den- 
sity at midrapidity by that obtained in inelastic p{p) + p 
collisions at the same collision energy, shown in Fig. O 
The similarity between the two data sets is remarkable. 
The inset in Fig. |31 shows the data with an expanded 
y-range and additional curved lines illustrating the ex- 
pectation for pure binary collision {Ncou) scaling at 200 
and 19.6 GeV. The dashed horizontal line represents the 
expectation for pure participant (Npart) scaling. In this 
representation, it becomes clear that the data at both 
energies follow a more iVpart-like dependence. 

Systematic errors dominate the charged particle den- 
sity measurements at 200 and 19.6 GeV. We find, how- 
ever, that most of these cancel in the ratio, leaving a base- 
line 3.0% overall uncertainty. This occurs as both analy- 
ses were performed with exactly the same method, detec- 
tor and with carefully matched centrality determinations. 
The main uncertainty comes from statistics of data and 
MC simulations (2.2%) and systematics in the primary 
charged particle detection efficiency (2.0%). Smaller sys- 
tematic contributions arise from background subtraction 
(0.4%) and uncertainty in the nucleon-nucleon cross sec- 
tion, (TjvAr (0.4%). We also include an additional cen- 
trality dependent systematic uncertainty that is largest 
in the more peripheral region and becomes negligible for 
the most central. This term takes into account the possi- 
bility that the estimated overall efficiency error may not 
entirely cancel in the ratio. 

As a final cross-check, the ratio of data at 200 to 19.6 
GeV, R200/19.6 is calculated in two distinct ways. First, 




FIG. 3: Centrality dependence of the measured Au + Au 
pseudorapidity density per participant pair divided by the 
corresponding value obtained in p{p) +p collisions. The light 
and dark error ellipses on the data points represent the sys- 
tematic errors associated with the measurement at 200 and 
19.6 GeV, respectively. There is an additional scale error as- 
sociated with the error on the value of the p{p) +p data points 
of 3.5% and 10% for the 200 and 19.6 GeV data, respectively. 
The figure inset has the same axes labels as the figure itself. 
The curved lines in the inset panel represent the binary colli- 
sion (Ncoii) scaling limit for both energies and the horizontal 
dashed line corresponds to pure Npart scaling. As in Fig. |5| 
error ellipses represent 90% C.L. limits. 



a more model-independent ratio was formed by dividing 
the data at each corresponding fraction of total interac- 
tion cross section. This is given by the solid squares in 
Fig. 21 (Au-fAu 1). Matching the centrality percentile 
bin at each energy, however, means there will be a differ- 
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FIG. 4: The ratio, R200/19.6, of the midrapidity pseudorapid- 
ity density per participant pair at 200 and 19.6 GeV versus 
(Npart)- The data binned by fraction of cross section are 
shown as closed squares and that binned by matching (Npart) 
are given as open squares. Also shown is the division of the 
inelastic p{p) +p collision data (open diamond) at Npart ~ 2. 
Curves give various calculations. The vertical error bars are 
combined statistical and systematic 1-a uncertainties. 



ence in the deduced (Npart) value at 19.6 and 200 GeV 
(see Table nj. In this case, the assigned (Npart) for each 
percentile bin given in the figure is taken as the average 
of the two individual (Npart) values. Second, the ratio 
was formed using a new set of centrality cuts for which 
each centrality bin width was varied, in an iterative fash- 
ion, in order to obtain bins at both 19.6 and 200 GeV 
that yield the same calculated average Npart- The data 
at both energies were then completely re-analyzed using 
this second set of centrality cuts. This result is given by 
the open squares of Fig. ^ (Au-|-Au 2). 

The results from the two types of ratio calculations 
are shown in Fig. 01 together with the predictions of 
two models and the two-component fit from Fig. [2 We 
find that both sets of data (closed and open squares) 
are in agreement, even within the significantly reduced 
systematic errors. This level of agreement gives fur- 
ther confidence that the systematic errors are cancel- 
ing as expected. Additionally, we find that the slope, 
and hence the centrality dependence, of both ratios is 
zero, within error. The most probable mean value of 
the (Au+Au 1) ratio data is found to be R200/19.6 — 
2.03±0.02(stat)±0.05(syst). We remind the reader that 
the ratio of 200 to 130 GeV data was found to be R200/130 
= 1.14±0.01(stat)±0.05(syst) 0. 

With the reduced systematic errors on the ratio now 
available, we return to a more detailed comparison of our 
results to calculations. As shown in Fig. 0] model calcu- 
lations predict quite different centrality dependences of 



R200/19.6 over the collision energy range of 19.6 GeV to 
200 GeV. We find that the Hijing calculation gives the 
expected increase of pQCD minijet production with cen- 
trality over this energy range, but the predicted increase 
is now in strong contradiction to the data. The fiat cen- 
trality dependence of the ratio is relatively well described 
by the parton saturation model calculation. 

In summary, PHOBOS has measured the charged par- 
ticle pseudorapidity density at midrapidity {\ri\ < 1) for 
Au + Au collisions at energies of ^^snn = 19.6 and 
200 GeV. We find an increase in particle production 
per participant pair for Au + Au collisions compared 
to the corresponding inelastic p{p) + p values for both 
energies. The ratio of the measured yields at 200 and 
19.6 GeV shows a clear geometry scaling over the cen- 
tral 40% inelastic cross section and averages to R200/19.6 
= 2.03±0.02(stat)±0.05(syst). A large increase in yield 
from hard processes, which contribute to multiplicity, is 
not apparent in the data, even over an order of magnitude 
range of collision energy. 
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